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■ Abstract This paper deals with the computation of the eigenvalues of two-parameter Sturm- 

Liouville (SL) problems using the Regularized Sampling Method, a method which has been effective 
in computing the eigenvalues of broad classes of SL problems (Singular, Non-Self- Adjoint, Non- 
Local, Impulsive,...). We have shown, in this work that it can tackle two-parameter SL problems 
with equal ease. An example was provided to illustrate the effectiveness of the method. 
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1 Introduction 

. In an interesting paper published in 1963, F. M. Arscott [2] showed that the method of separation of 

variables used in solving boundary value problems for Laplace's equation leads to a two-parameter 
eigenvalue problem for ordinary differential equations with the auxiliary requirement that the so- 
lutions satisfy boundary conditions at several points. This has led to an extensive development of 
multiparameter spectral theory for linear operators (see for instance [3-4], [7-10], [12], [15], [30-33], 
[35], [37], [39-42]). In the paper [12], the authors give an overview results on two-parameter eigen- 
value problems for second order linear differential equations. Several properties of corresponding 
1 eigencurves are given. In [15] the authors have obtained interesting geometric properties of the 

eigencurves (for instance transversal intersections is equivalent to simplicity of the eigenvalues in 
the sense of Chow and Hale). All the above works are concerned with the theoretical aspect of 
existence of eigenvalues. Also, several authors have dealt with the theoretical numerical analysis 
of two-parameter eigenvalue problems (see [5], [11], [13], [34], [36], [38] and the references therein). 
Eigenvalue problems have played a major role in the applied sciences. Consequently, the problem of 



oo 

' computing eigenvalues of one-parameter problems has attracted many researchers (see for example 

[1, 6, 16,17,19-29] and the references therein). 

Concerning the computations of eigenvalues of one-parameter Sturm-Liouville problems, the 
authors in [16] introduced a new method based on Shannon's sampling theory. It uses the analytic 
properties of the boundary function. The method has been generalized to a class of singular 
problems of Bessel type [17], to more general boundary conditions, separable [19] and coupled [21], 
to random Sturm-Liouville problems [23] and to fourth order regular Sturm-Liouville problems 
[28]. The books by Atkinson [4], Chow and Hale [30] , Faierman [33] , McGhee and Picard [37], 
Sleeman [39] , Volkmer [42] and the long awaited mongraph by Atkinson and Mingarelli [5 contains 
several results on eigenvalues of multiparameter Sturm-Liouville problems and the corresponding 
bifurcation problems. However, no attempt has been made to compute the eigenvalues of two- 
parameter Sturm-Liouville problems using the approach based on the Regularized Sampling Method 
introduced recently by the first author in |25j to compute the eigenvalues of general Sturm-Liouville 
problems and extended to the case of Singular [24 , Non-Self- Adjoint [23], Non-Local [20], Impulsive 
SLPs [22]. [21]. We shall consider, in this paper the computation of the eigenpairs of two-parameter 
Sturm-Liouville problems with three-point boundary conditions using the Regularized Sampling 
Method. 

2 The Characteristic Function 

Consider the two-parameter Sturm-Liouville problem 
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/ -y" + qy = (nlm + i£w 2 )y , o < x < 1 , , 

\ y(0) = , y(c) = , »(1) = 

where w\, w 2 are positive and in C 2 [0, 1] and g € L[0, 1] and c G (0, 1) some given constant. 

By an eigenvalue of (|2.1[) we mean a value of the couple for which problem (|2.1[) has 

a non trivial solution. Conditions that insure the existence of eigenvalue arc given in [3], [9], [10], 
[15], [30], [31], [40] and [41]. In fact, under fairly general conditions it has been shown (see [15], [30]) 
that there are smooth curves of eigenvalues (actually eigenpairs). Our objective is to effectively 
localize the eigencurves in the parameter (pi, fi 2 )-p\ane.We should point out that we have restricted 
our attention to Dirichlet boundary conditions in order to eliminate technical details that might 
obscur the ideas. 

We shall associate to (|2.1[) the initial value problem 

/ y" + (Mi^i + ^\w 2 )y = qy , < x < 1 , . 

\ y(0) = , j/(0) = 1 [Z - Z) 

and deal first with the unperturbed case (q = 0) then with the perturbed case (q 0). 



2.1 The unperturbed case (q = 0) 

In this case, (|2.2[) reduces to 

ip'{ + wtpi = , < x < 1 

¥>i(0) = 0, <^i(o) = i 

where w = ^\w\ + fJ, 2 w 2 . 



(2.3) 



Theorem 1 The solution ipi of 12. 3\) is an entire function of (fix , p 2 ) £ C 2 for each fixed x € (0, 1] 
of order (1, 1) and type (ai(x), <J 2 {x)) and satisfies the estimate, 

\f\{x)\ < Kiexp[ai(x) |/xi| + a 2 (x) |/x 2 |] , (MijMa) € C 2 

/or eac/i fixed x 6 (0, 1] where, o~i{x) — 2 {x J Q X Wi(£)<i£ } 2 /or £ = 1,2. 
Proof. From (|2.3[) we get the integral equation 

Vi{x)=x- I (x - [rfMO + »lM0] MZH (2-4) 

(2.5) 

(2.6) 



Let 

(/3i,„+i(a;) = $*{x - [f4u>i(0 + (4w 2 {£)] (pi, n (0<^ ,n>0 
We shall show, by induction on n, that 



\(Pi,n(%)\ < 



n!(n + l)! t Jo 



|/*i|V(0 + M 2 «*(0 d£\ ,n>0 



It is true for n = 0. Assume it is true for n. We shall show that it is true for n + 1. Indeed, from 
(12.51). we have 



|¥?i,n+i(aOI < / (s-C) iMil ^i(0 + |M2| u> 2 (£) 



n!(n + l)! Py o 



|Mi| 2 Wi(t) + |^2| 2 W2(t) dr ^ d£ 



(2.7) 
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Using the fact that the expression (x — £)£ attains its maximum at £ = ^z^x, we get 



| <£>i,„+i (a;) | < 



1 1 



< 



n!(n + 1)! n + 1 n + 2 \n + 2 

x 



(n- 




r r r 















|//i| 2 wi(r) + |« 2 | 2 w 2 (r) 



n+l 



(n+l)!(n + 2)! 



Ml| 2 WlM + |/i2| 2 W2(t) 



n+l 



(2.8) 



that is (|2.6[) is true for n+l. Hence, it is true for all n > 0. 

Now, </?i(a:) = X)n>o( — l) n yi,ra( a; ) ano - the series is absolutely and uniformly convergent since 



l¥>i(z)l 



^(-l)>i,„( a 

n>0 



n>0 



n>0 



n!(n + 1)! 



|ui| 2 u;i(r) + \H2\ 2 w 2 (t) 



(It 



xhU2Jx \^ 1 \ 2 w 1 {t) + \^ 2 \ 2 w 2 {t) dr})x 



(It 



(i 



\[ii\ 2 wi(r) + u 2 | 2 w 2 (r) 



where ii is the modified Bessel function of the first kind order 1 

■ z \ 2n+l 

nUn 4- 11! V 

n>0 



/l(z) ~^n!(n + l)! V 2 



Using the fact that Ji(z) 



/27TZ 



as z — > 00, we get 
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\<pi(x)\ < Kiexp 

< i^iexp I +ct 2 (x) |u 2 |] 



|«l| 2 + |M2| 2 W 2 (t) 



(It 



(2.9) 



(2.10) 



Therefore, tp\ is an entire function of (/ii,u 2 ) G C 2 , as a uniformly convergent series of entire 
functions, for each fixed x S (0, 1], of order (1, 1) and type (oi(x), cr 2 (a;)).This concludes the proof. 



We shall make use of the Liouville-Green's transformation 



to bring (|2 .2[) to the form 



z(t) = {w(x)}* ipi{x) 



f |f + {1 + r(<)}z = , < i < Ji y/Mg)d£ 

\ «(o) = o l |(o) = Wo)}-i 

which can be written as an integral equation 

z(t) = {w(0)}~^ sini - f sin(t - T)r(T)z(T)dr 
Jo 

where r(t) = {w(x)} 4 £j 4 

' J \x—x(t) 



(2.11) 



(2.12) 



(2.13) 
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Returning to the original variables, we deduce that (fi satisfies the integral equation 
^(x) = {w(0)w(x)y^ sin j^ - J sin (J Vw(O d O*P(x, x)<Pi (x)dx (2.14) 



where 



ip(x, x) = {w(x)} 4 {w(x)Y 



\w"{x) {w(x)} ' + 4 HOe)} 2 {w^)}" 



(2.15) 



4 16 
We shall present next some estimates whose proofs are immediate and left to the reader. 
Lemma 2 The function i[)(x,x) satisfies the estimate 

\ip(x,x)\ ~ - — — - as \m\ + \fj, 2 \ -> oo, (mi,A*2) G R 2 (2.16) 
iMil + |M2| 

Lemma 3 TTie function a defined by 

a(x) — (sine {cr 1 (x) fix + cr 2 (x)/U 2 }) m 

where sinc(z) = z sinz and m is a positive integer, is an entire function of (/Ui,/tx 2 ) £ C 2 /or 
eac/i fixed x G (0, 1] of order (1, 1) and £?/pe (<7i(x), <t 2 (x)). Furthermore, a satisfies the estimate 

H 1 )! ~ i — rn — i' as ImiI + Ims I ->■ oo, (mi,m 2 ) e k 2 . 

iMil + IM2I 

Lemma 4 TTie function ipi satisfies the estimate 

M 1 )! ~ i M 4 , i > as M + M ^00,(^,^3) el 2 (2.17) 

\m\ + |M2| 

Combining the above results, we obtain the following theorem, 

Theorem 5 The function oup\ is an entire function of (fii,fi 2 ) £ C 2 for each fixed x e (0, 1] of 
order (1,1) and type ((to + l)cr 1 (x), (m + l)er 2 (x)) and satisfies the estimate 



K(x) 

(m + m) ? 



|a(x)^i(x)| - — ; — — as ^ + |^2 1 -> 00, (/ii,A*2) G 



where X depends on x e (0, 1] but is independent of (/xi,/x 2 ) • 
Let PW/3 lt /3 2 denote the Paley- Wiener space, 

/ h(zi, z 2 ) entire / |ft(zi,2 2 )| < Cexp{/3i |zi| + /3 2 |z 2 |} , \ 

we have, 

Theorem 6 a(x)<^i(x) , as a function of (/xi,/x 2 ) belongs to the Paley-Wiener space PW^ 1 ^ 2 
where (f3\, (3 2 ) = ((m + l)eri(x), (to + 1)ct 2 (x)) for each fixed x E (0, 1]. 
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2.2 The perturbed case (q 7^ 0) 



Let if 1 , <f 2 be two linearly independent solutions of f" + wip — satisfying ipi(0) = ^(O) = , 
ip'i(0) = y2(0) = 1 then the method of variation of parameters shows that (|2.2I) can be written as 
the integral equation 



y(x) = ip 1 (x)+ I Hx^MtmW ( 2 - 18 ) 
where ®(x,g) = ifi(£)<f 2 (x) - ^(OVifa)- 







Here again, it is not hard to show that y(x) is an entire function of (/zi,/i2) for each x G (0, 1], 
of order (1, 1) and type (<7i(x), 02(2;)). Multiplication by a gives a function a(x)y{x) of (/^i,/^) in 
a Paley- Wiener space PWp 1 ^p 2 for each x G (0, 1]. More specifically, we have the following, 

Theorem 7 The function y defined by y(x) — a(x)y(x) , belongs to PW^r 2 where (/3 X , /3 2 ) = 
((to + l)<7i(x), (to + l)a~2(x)) as a function of (/ii, /i 2 ) G C /or eac/i x G (0, 1], and satisfies the 
estimate, 

K(x) 

\a{x)y(x)\ ~ - — ■ ■ — -t^+t, as + |/i 2 | -> 00, (/xi,/i 2 ) G K 2 . 

(M + \m) 

where K depends on x G (0, 1] 6wi is independent of (pi, (12) ■ 

Proof. Since $ ra + w$ = , <&(t, t) = and 3>x(t, i) = 1, we have, 
$(x,t) = {w(t)w(2;)}~* sin (J y/w(£)dtj 

i^J V w (0dij i>(x,x)$(x,t)dx (2.19) 



sm 



so that, 

Mx,t)\ ~ , ,^ 4 , ; <K 5 , as |Mi| + |M2| oo,(mi,M2) eK 2 . (2.20) 

M + IM2I 

from which we get, after using Gronwall's lemma on (|2.18p and the estimate for fi, 
\y(x)\ < Ki exp [<ji(x) [/xx I + 02(2;) I/X2I] exp Ik 5 J \q{t)\di) 
12/(1)1 < tf 7 exp[<7i(l) +o- 2 (l) |// 2 |] , as |^i| + | M2 | -»• 00, (jjh, M2) G K 2 . (2.21) 

and 

Hl)y(l)| < A' 8 exp [(to + l)o-i(l) |/ji| + (m + l)o- 2 (l) |/J 2 |] , as \ni\ + \/j, 2 \ ->■ 00, (/xi,£t 2 ) G R 2 . 

(2.22) 

Furthermore, we have, 

if (a;) 

|a(x)y(a;)| ~ - — ■ ■ — -^+r,as \fii\ + \fi 2 \ -> 00, (^1,^2) G M 2 . 

(M + IM2I) 

where K depends on x G (0, 1] but is independent of (jJ,i,fjL 2 ) ■ ■ 

To summarize, in both cases, unperturbed and perturbed, the transform y(x; fii, fx 2 ) of the solu- 
tion y(x; hi, (12) of (I2.2p is in a Paley- Wiener space PWp 1 ^p 2 where (/?i, ^2) = ((to + I)cti(x), (to + 
l)<72(x)). Thus y(x\ /Hi, ^2) can be recovered at each x G (0,1] from its samples at the latice 
points (Mlj,M2fc) = (j (m+ i] ffl ( I ) > fc (m+iWz) )> e Z2 usin S the rectangular cardinal series 

(0315131511), 
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Theorem 8 Let f e Ff ftA tften 

,/„ \ v ,/„ „ N sin/3i(/xi - m,) sin ^ 2 (^ 2 - ^ 2fc ) 

f/ie convergence of the series being uniform and in L^ id/ ^ 2 (M. 2 ), and \i mn = nir/(3 m , m = 1,2, 

nez. 



Let Cll = Cl(l), (7 2 1 = ^(l), CTl2 = (71 (c), (7 22 = CT 2 (c). 

The eigenpairs are therefore (/z^/xl) where (^1,^2) solve the nonlinear system 



?S M ! = n (2.23) 
S 2 (/ii,M2) = 



where, 



{D ,, \ A 1 V^°° ~n . ,, ,, N sin2(m+l)g-ii(/Ji-^ij) sin 2(m+l)g 2 i (p2 -M2fc ) 

£>l^l,/!2j - Q (i) 2^j=-oo 2^fe=-ooi'l 1 'My I A i 2fcJ 2(m+l)<7ii(Mi-A»ii) 2(m+l)a 2 i (ti2~H2k) 
R f/, i< ^ A _1 v^oo ^p 00 JTfV ,/ N sin2(m+l)ji 2 (Mi-Mij) sin 2(m+l)g 22 (p 2 -^ 2fc ) 

£»HMl,A*2j - ^Z^j=-cx>Z^fc=-cx>yV C >Ml?,M2fcJ 2(m+l)CTi 2 (fti-;*ij) 2(m+l) C r 22 ( M2 - Al2fc ) 

(2.24) 

3 A numerical example 



We shall consider in this section the two-parameter Sturm-Liouville problem with three- point 
boundary conditions given by 

-y" = + ii\x)y , < x < 1 , , 

2/(0) = 2/(0.7) = y(l) 

The general solution y of the first differential equation can be expressed in terms of Ai and Bi 
functions and their first dirivatives as 

(-DV. (a, Bi (i3|pM) _ Ai (MM) Bi (JEM 

y(x;fj, 1 ,fj, 2 ) = 7 7 x 7 x 7 x 7 

(3.2) 

Thus the eigenpairs (/xf, jLt|) can be obtained from the solutions (111,112) of the system 

2/(1; Mi, M2) = , 33 ^ 
2/(c;/ii,/i 2 ) =0 

For numerical purposes we have truncated the associated scries to \k\ < N — 50 and took 
m = 5 in the function a. Thus, the approximate eigenpairs are seen as solutions of the system 



27^(1; Mi, M2) = 
yW(c; Ml ,M 2 )=0 



[N}(„... ,^_n ( 3 - 4 ) 



The next table shows the exact eigenpairs together with their approximations using the Regularized 
Sampling Method (RSM). 
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ji\ (exact) 


jj,2 (exact) 


Mi (RSM) 


[i2 (RSM) 


7.788149097670813 


7.5908224365786845 


7.78814883048352 


7.590823605399021 


4.194056047341936 


22.273796542861913 


4.194057357011064 


22.273795699572364 


15.597607295800684 


15.165889997583099 


15.597605904429917 


15.1658929384375 


13.87615754699624 


30.597837626955204 


13.876157718920586 


30.59783713797321 


23.40242244972004 


22.744341450749697 


23.40242640980619 


22.744329536815076 


9.51130982713563 


44.296491367611544 


9.511333341150682 


44.29647446214639 


39.31833491621106 


15.670248767892955 


39.31833152483279 


15.670257662031721 


47.18632817417246 


24.620489734469363 


47.186345272524804 


24.62038504864854 


46.81208951678993 


45.483251580157955 


46.813097550330646 


45.480149444567346 


30.03437464767369 


46.558148101315844 


30.034392317076748 


46.55811263809656 



4 Conclusion 

In this paper, we have successfully computed the eigenpairs of two-parameter Sturm-Liouville prob- 
lems using the regularized sampling method. A method which has been very efficient in computing 
the eigenvalues of broad classes of Sturm-Liouville problems (Singular, Non-Self- Adjoint, Non- 
Local, Impulsive,...). We have shown, in this work that it can tackle two-parameter SL problems 
with equal ease. An example was provided to illustrate the effectiveness of the method. 

Acknowledgements The authors are grateful to KFUPM for its usual support. 
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